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Abstract. New estimates of the distances of 36 nearby 
galaxies are presented based on accurate distances of 
galactic Cepheids obtained by Gieren, Fouque and Gomez 
(1998) from the geometrical Barnes-Evans method. 

The concept of 'sosie' is applied to extend the distance 
determination to extragalactic Cepheids without assum- 
ing the linearity of the PL relation. Doing so, the distance 
moduli are obtained in a straightforward way. 

The correction for extinction is made using two photo- 
metric bands {V and /) according to the principles intro- 
duced by Freedman and Madore (1990). Finally, the sta- 
tistical bias due to the incompleteness of the sample is cor- 
rected according to the precepts introduced by Teerikorpi 
(1987) without introducing any free parameters (except 
the distance modulus itself in an iterative scheme) . 

The final distance moduli depend on the adopted ex- 
tinction ratio Rv / Ri and on the limiting apparent magni- 
tude of the sample. A comparison with the distance mod- 
uli recently published by the Hubble Space Telescope Key 
Project (HSTKP) team reveals a fair agreement when the 
same ratio Ry/Ri is used but shows a small discrepancy 
at large distance. 

In order to bypass the uncertainty due to the metallic- 
ity effect it is suggested to consider only galaxies having 
nearly the same metallicity as the calibrating Cepheids 
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(i.e. Solar metallicity). The internal uncertainty of the dis- 
tances is about 0.1 magnitude but the total uncertainty 
may reach 0.3 magnitude. 

Key words: galaxies: distances and redshift - galaxies: 
stellar content - cosmology: distance scale 



1. Introduction: Discussion of the problems 
related to Cepheids 

As an extension of our study of the kinematics of the lo- 
cal universe (KLUN-I-) we need an accurate value for the 
global Hubble constant and accurate distances of individ- 
ual galaxies. The calibration of the distance scale is thus 
a fundamental step in this process. The aim of this work 
was to calibrate the distance scale from nearby galactic 
Cepheids for which the HIPPARCOS satellite measured 
geometrical parallaxes. This should avoid the step of cal- 
ibrating the distance scale by assuming a given distance 
to the Large Magellanic Cloud (LMC). Unfortunatelly, it 
turns out that these measurements are very difficult to 
use due to a statistical bias (Lutz and Kelker, 1973). The 
difficulties can be solved by proper treatment, like the 
one proposed by Feast and Catchpole (1997). It has been 
shown that this leads to unbiased results (Pont et al., 1997; 
Lanoix et al. 1999), 

On the other hand, individual measurements of 
Cepheids from HIPPARCOS are relatively inaccurate be- 
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Table 1. Slopes of the PL relation. 



source 


av 


ai 


GFG(MW) 


-3.037 ±0.138 


-3.329 ±0.132 


GFG(LMC) 


-2.769 ± 0.073 


-3.041 ± 0.054 


OGLE(LMC) 


-2.765 


-2.963 



cause of the distance of galactic Cepheids. Excluding a 
UMi which does not pulsate in the fundamental mode, 
the best geometrical parallax of an individual Cepheid 
obtained from HIPPARCOS is 3.32± 0.58 marcsec for 
6 Cephee. This leads to an uncertainty in the distance 
modulus of 0.38 magnitude. In comparison, the quasi- 
geometrical method of Barnes-Evans applied to Cepheids 
(Gieren, Fouque, Gomez, 1998; hereafter GFG), gives dis- 
tance moduli with a typical uncertainty less than 0.1 mag- 
nitude (the external error can be estimated to about 0.2 
magnitude according to Table 7 in GFG). We call this 
method quasi-geometrical because it requires only a few 
assumptions. The method is independent of any determi- 
nation of the LMC distance and has a relatively small sys- 
tematic error (about 0.2 magnitude). Thus, we decided to 
calibrate the distance scale using the work done by Gieren, 
Fouque and Gomez (1998). 

Nevertheless, other difficulties appear. The slope of the 
Period-Luminosity relation (hereafter, PL relation) deter- 
mined from the adopted calibrating galactic Cepheids dif- 
fers from the slope obtained for the LMC by the same 
authors (GFG) (Table 0).. For the LMC, the slopes in 
V and I bands are now confirmed by the OGLE survey 
(Udalski. et al., 1999). What slope should we adopt? 

The true physical relation is actually a Period- 
Luminosity-Color (hereafter, PLC) relation written as 
M = alogP ± (iCo ± 7, where M is the absolute mag- 
nitude and Co the intrinsic color. The PL relation is sim- 
ply the projection of the PLC onto the P-L plane. In the 
PLC relation the slope dM/dlogP is constant. However, 
the observed slope of the PL relation depends on the dis- 
tribution of observed Cepheids in the PLC plane (i.e., on 
the color distribution of the sample). Hence, the slope in 
a given photometric band may partially depend on the 
metallicity, because it affects the intrinsic color. Linear 
non-adiabatic models do predict that the slope is con- 
stant when one uses bolometric magnitudes (Baraffe et 
al., private communication), whereas non- linear models 
predict that the slope depends on the metallicity also for 
the bolometric magnitudes (Bono et al., 2000 and refer- 
ences therein) and predict that the slope in a given band 
depends on the metallicity. Because the metallicity of the 
LMC differs from the metallicity in the Solar neighbour- 
hood, the choice of slopes in different bands is difficult. In 
order to avoid this dilemma we decided to apply the method 



of 'sosie' (Paturel, 1984) because it does not require knowl- 
edge of the slope and zero point of the PL relation Q 

The correction for extinction produced by interstellar 
matter is another difficulty. It can be solved by assuming 
that the extinction law is universal. We will thus assume 
that the extinction on an apparent magnitude is propor- 
tional to the color excess {A\ = R\{C — Co), where C is 
the reddened color). The factor of proportionality Rx is 
taken from tabulations (e.g., Cardelli, Clayton & Mathis, 
1989 ; Caldwell & Coulson, 1987 ; Laney & Stobie, 1993). 
It depends on both the considered band and color. With 
such an assumption it is possible to use the Freedman 
and Madore (1990) precepts of de-reddening. Two bands 
are needed in order to calculate a color. Because most ex- 
tragalactic Cepheids are measured in V- and I-band from 
The Hubble Space Telescope (hereafter, HST), we will use 
these two bands. Thus, the Freedman and Madore (1990) 
de-reddening method will be adapted to the sosie method, 
used in V and I photometric bands. 

Finally, an ultimate difficulty comes from the incom- 
pleteness bias. This bias was first studied by Teerikorpi 
(1987) for application to galaxy clusters (Bottinelli et al., 
1987). It was ffist denounced by Sandage (1988) in appli- 
cation to the PL relation and re-discussed later by Lanoix, 
Paturel and Garnier (1999a). The sample to which we are 
applying the PL relation must be statistically represen- 
tative of the calibrators themselves. Indeed, due to the 
intrinsic scatter of the PL relation, there is a given dis- 
tribution of absolute magnitudes at a given period. At 
increasing distances the fainter end of this distribution 
is progressively missed and the distribution of the actual 
sample changes. Restricting the sample to Cepheids with 
a period larger than a given limiting period reduces this 
bias. The limiting period depends on a first estimate of 
the distance, on the apparent limiting magnitude and on 
the characteristics of the PL relation (dispersion, slope 
and zero-point). In fact, the full theory of Teerikorpi is 
applicable. The method is much more complete than the 
rough rule of thumb used as a quick approach in an ap- 
plication in which a detailed treatment was not needed. 
However, we want to derive final distance moduli and the 
precise bias correction must be used. Note that the slope 
and zero point of the PL relation are needed but only as 
second order terms and thus, the uncertainties mentioned 
about their choice do not present any significant difficulty 
(this will be confirmed in section 4.3). The incompleteness 
bias will be corrected using the precepts given by Teerikorpi 
(1987). 

In section 2 we will describe the material used for this 
study: the calibrating sample by GFG and our extragalac- 
tic Cepheid database (Lanoix et al., 1999b). 



this method was first introduced to solve the same kind of 
problems for the TuUy-Fisher relation (1977) 
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In section 3 we describe the 'sosie' method and give the 
basic equation for the calculation of the distance modulus 
of an extragalactic Cepheid. 

In section 4 we give the results obtained for 1840 
Cepheids belonging to 36 nearby galaxies described in the 
previous section. We also discuss these results and com- 
pare them with those recently published by Freedman et 
al. (2001). 

2. Observational material 

The guideline in the constitution of the observational ma- 
terial is the selection of the most secure observations. This 
leads us to reject some data, as explained below, both 
galactic and extragalactic. 

2.1. The list of galactic Cepheids 

The starting point of our study is the choice of the galac- 
tic Cepheids used for the calibration. We adopt the list 
given in Gieren, Fouque and Gomez (Table 3 in GFG) 
but we rejected three Cepheids (EV Set, SZ Tau and 
QZ Nor) because they do not pulsate in the fundamen- 
tal mode (they are overtone Cepheids). They correspond 
to the three lowest periods of the list. Because we use 
only the V and I photometric bands, three Cepheids are 
also rejected (CS Vel, GY Sge and S Vul) because they do 
not have I-band magnitude. Thus, 28 Cepheids remain. 
Their distance moduli are adopted directly from Table 5 
given by GFG. Only three Cepheids have a mean error 
in their distance modulus larger than 0.1 magnitude. We 
give in Table ^ the adopted calibrating sample of galactic 
Cepheids. 

2. 2. The list of extragalactic Cepheids 

In 1999 we have constructed an Extragalactic Cepheid 
database (Lanoix et al., 1999b) by collecting 3031 pho- 
tometric measurements of 1061 Cepheids located in 33 
galaxies. This list has been updated. Especially, the V 
and I band measurements by Udalski et al. (OGLE sur- 
vey, 1999) were added for the LMC from the data avail- 
able through 'astro/ph9908317'. The new database con- 
tains 6685 measurements for 2449 Cepheids in 46 galax- 
ies. In order to make this compilation available, the full 
contents of the extragalactic part will be published in elec- 
tronic form for the A&A archives at CDS. A description 
is given in the Annex. 

In this database, each light curve has been inspected in 
order to describe the main features. In the present study 
only light curves considered as 'Normal' are used |^ We 

^ Lanoix et al. give eight classes of light curves: 'Normal', 'Sy- 
metrical', 'Bumpy', 'Scattered', 'Overtone', 'Low amplitude', 
'Peculiar', 'No curve'. A 'Normal' light curve is characterized 
by a non-symetrical variation: a fast increase and a slower de- 
crease. 



Table 2. Adopted calibrating sample of galactic Cepheids. 
Column 1: Name of the galactic Cepheid; Column 2: 
log of the period (P in days); Column 3: adopted dis- 
tance modulus and its mean error according to Gieren, 
Fouque and Gomez (1998); Column 4: Mean V-band 
apparent magnitude; Column 5: Mean I-band apparent 
magnitude. 



Cepheid 


logP 


fi ± m.e. 






BF Oph 


0.609 


9.50± 


0.11 


7.33 


6.41 


T Vel 


0.666 


10.09± 


0.02 


8.03 


7.01 


CV Mon 


0.731 


10.90± 


0.05 


10.31 


8.68 


V Cen 


0.740 


9.30± 


0.02 


6.82 


5.81 


BB Sgr 


0.822 


9.24± 


0.02 


6.93 


5.84 


U Sgr 


0.829 


8.87± 


0.01 


6.68 


5.45 


S Nor 


0.989 


9.92± 


0.03 


6.43 


5.41 


XX Cen 


1.039 


10.85± 


0.06 


7.82 


6.75 


V340 Nor 


1.053 


11.50± 


0.13 


8.38 


7.15 


UU Mus 


1.066 


12.26± 


0.09 


9.78 


8.49 


U Nor 


1.102 


iU. M± 


U.O 1 


9.23 


7.36 


JjlN r^up 


l.loD 


12.92± 


0.05 


O GO 


O.OO 


LS Pup 


1.151 


13.73± 


0.04 


10.45 


9.06 


VW Cen 


1.177 


13.01± 


0.04 


10.24 


8.77 


VY Car 


1.277 


11.42± 


0.04 


7.46 


6.28 


RY Sco 


1.308 


10.47± 


0.04 


8.02 


6.30 


RZ Vel 


1.310 


11.17± 


0.03 


7.09 


5.85 


WZ Sgr 


1.339 


11.26± 


0.02 


8.02 


6.53 


WZ Car 


1.362 


12.98± 


0.14 


9.26 


7.95 


VZ Pup 


1.365 


13.55± 


0.04 


9.63 


8.28 


SW Vel 


1.370 


11.99± 


0.06 


8.12 


6.83 


T Mon 


1.432 


10.58± 


0.07 


6.12 


4.98 


RY Vel 


1.449 


12.10± 


0.05 


8.37 


6.84 


AQ Pup 


1.479 


12.75± 


0.04 


8.67 


7.12 


KN Cen 


1.532 


12.91± 


0.06 


9.85 


7.99 


L Car 


1.551 


8.94± 


0.05 


3.73 


2.59 


U Car 


1.589 


11.07± 


0.04 


6.28 


5.05 


SV Vul 


1.654 


12.32± 


0.07 


7.24 


5.75 



reject all peculiar light curves including light curves clas- 
sified as 'low amplitude' because they are often associated 
with overtone Cepheids. 

Only the mean V and I band magnitudes are kept. 
When several magnitudes are averaged from different 
sources we keep the mean only if the mean error is less 
than 0.05 magnitude. It is to be noted that HST mea- 
surements of seven galaxies |^ have been analyzed by two 
independent groups. This leads to two different sets of 
magnitudes. Independent treatment of both sets shows 
that the distance modulus differs by less than 0.1 mag- 
nitude, except for IC4182 for which the difference is 0.28 
magnitude (Lanoix, private communication). Because we 
have no means to decide which set is the best we decided 
to keep them both. 

The final catalogue (Table ^ results in 1840 extra- 
galactic Cepheids. They belong to 36 galaxies, 27 of which 

^ IC4182, NGC3368, NGC3627, NGC4496A, NGC4536, 
NGC4639 and NGC5253 
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Table 3. Sample of extragalactic Cepheids. Column 
1: Name of the host galaxy; Column 2: Name of the 
Cepheid according to the following reference; Column 3: 
Reference (coded) from which the Cepheid name is taken; 
Column 4: log of the period (P in days); Column 5: 
Mean V-band apparent magnitude; Column 6: Mean I- 
band apparent magnitude. Only a part of the table is 
given. The rest is available in electronic form. 



galaxy 


Cepheid 


Ref. 


log? 


{V) 




IC4182 


Cll 


Gib99 


1.423 


23.10 


22.21 


LMC 


109838 


Uda99 


0.732 


16.14 


15.11 


NGC1365 


V32 


Sil98 


1.460 


26.77 


25.94 


NGC1425 


C15 


Mou99 


1.295 


26.63 


25.90 


NGC2090 


C13 


Phe98 


1.461 


25.44 


24.55 


NGC224 


FI13 


Fre90 


1.497 


19.24 


18.33 


NGC2541 


C25 


Fer98 


1.270 


25.68 


24.90 


NGC3031 


C13 


Fre94 


1.270 


23.56 


22.75 


NGC3109 


P2 


Mus98 


0.722 


22.18 


21.87 


NGC3198 


C19 


Kel99 


1.220 


26.23 


25.12 


NGC3319 


C13 


Sak99 


1.398 


25.61 


24.89 


NGC3351 


C25 


Gra97 


1.207 


25.77 


24.49 


NGC3368 


C09 


Gib99 


1.483 


25.13 


24.11 


NGC3621 


C14 


Raw97 


1.498 


23.28 


22.76 


NGC3627 


C14 


Gib99 


1.366 


24.66 


23.46 


NGC4258 


MA014 


Mao99 


1.330 


24.65 


23.88 


NGC4321 


C9 


Fer96 


1.700 


25.93 


24.88 


NGC4414 


CI 


Tur98 


1.658 


25.89 


24.85 


NGC4496 


C24 


Gib99 


1.717 


25.27 


24.26 


NGC4535 


C35 


Mac99 


1.390 


26.14 


25.22 


NGC4536 


C12 


Gib99 


1.484 


25.81 


24.89 


NGC4548 


C09 


Gra99 


1.270 


25.96 


25.38 


NGC4603 


2984 


New99 


1.570 


27.19 


26.37 


NGC4639 


C14 


Gib99 


1.717 


26.33 


25.28 


NGC4725 


C09 


Gib98 


1.590 


24.85 


23.87 


NGC5253 


C07 


Gib99 


1.025 


23.71 


22.86 


NGC5457 


V4 


Kel96 


1.471 


23.51 


22.78 


NGC598 


V31 


Chr87 


1.572 


19.17 


18.14 


NGC7331 


V4 


Hug98 


1.354 


26.13 


24.93 


NGC925 


V18 


Sil96 


1.439 


24.99 


23.97 


SEXB 


V2 


Sa85b 


1.444 


20.60 


20.00 



come from HST observations and 9 from ground-based ob- 
servations. The full Table is available in electronic form in 
the A&A archives at CDS. 



3. Method of sosie 

The method of 'sosie' was introduced (Paturel, 1984) to 
avoid the problem encountered in the practical use of the 
TuUy-Fisher relation (TuUy and Fisher, 1977), a linear re- 
lation between the absolute magnitude of a galaxy and its 
21-cm line width. Here we are in similar conditions with 
a linear relationship between the absolute magnitude and 
an observable parameter, the logarithm of the period. In 



French, the word 'sosie' refers to someone who looks very 
similar to someone else without being necessarily genet- 
ically related. Here two Cepheids will be considered as 
'sosie' if their light curves have the same shape and if 
they have the same period (within a given error) . Because 
of the selection based on the shape of the light curve we 
will consider that all Cepheids of our sample pulsate in the 
fundamental mode. They all obey the same P-L relation. 

We write the distance modulus of a calibrating 
Cepheid and of an extragalactic Cepheid through a uni- 
versal PL relation. The calibrating Cepheid is identified 
with subscripts 'c' and no subscript for the extragalactic 
one. Presently, we assume that both stars have the same 
metallicity and the same intrinsic color. We will see how 
to bypass this problem, later. 



= m° - a log Pc - b, 
fi — ni° — a log P — b 



(1) 

(2) 



rn° is the apparent mean magnitude in a given band. The 
superscript 'o' means 'corrected for extinction'. If one se- 
lects an extragalactic Cepheid having the same period as 
the calibrating one, i.e., logP = log Pc, the distance mod- 
ulus of the extragalactic Cepheid is then 



(3) 



The distance modulus of the extragalactic Cepheid is de- 
duced without having to know the slope and zero-point of 
the PL relation. 

In order to correct for extinction we apply the previous 
equation to two different bands and express the extinction 
term as a function of the color excess E = E{B — V). In 
order to make the notations clearer we note the apparent 
magnitudes V and / for the two considered bands. From 
Eq. ^ one has: 



H = fic + V- RvE -Vc + RvEc 
/i = /ic + / - RiE + RiEc 
which can be written as 
fi = fic + V-Vc-Rv{E-E,) 
^l = ^lc + I-Ic-RI{E-E,) 



(4) 
(5) 

(6) 
(7) 



Then, eliminating E — Ec between the two previous equa- 
tions we obtain: 



1 



Ri 



(8) 



This is the desired equation. It can be written in a more el- 
egant manner by using the reddening-free Wesenheit func- 
tion (Van den Bergh, 1968): 



W = 



V^{Rv/Ri).I 
1 - {Rv/Ri) 



(9) 
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^ = ^^^ + W -Wc 



(10) 



In practice, the intrinsic color is not known and this equa- 
tion is valid only for Cepheids of the same intrinsic color 
and metallicity. Thus, for a true sample, we will write it 
as (see the discussion below): 



(m) = Mc + {W- W,) 



(11) 



W is an observable quantity. Then, the mean distance 
modulus of a sample of Cepheids which have the same pe- 
riod of pulsation as a calibrating Cepheid can be obtained 
directly from Eq. |ll|. 

The physical relationship in this result is a Period- 
Luminosity-Color relation. This means that we should 
search for sosie of calibrators by considering both their 
similarity in logP and intrinsic color Co- But the intrinsic 
color is not observable. Thus, equation |ll| must be consid- 
ered as a statistical relation exactly as the PL relation. Be- 
cause of the statistical relation between Co and logP, the 
selection in log P will guarantee that a calibrator Cepheid 
and its sosies have, on average, the same intrinsic color. 
So, the problem of the intrinsic color is partially bypassed. 
For the metallicity problem, the solution is to consider 
that the method is valid only for galaxies having nearly 
the same metallicity as the calibrating Cepheids. In the 
present paper this means that, stricto sensu, only galaxies 
with a nearly Solar metallicity can be considered as valid. 
In practice, we applied the method to different kinds of 
galaxies without noting strong metallicity dependence. 

As a test, we apply the method to the calibrating 
sample itself. Indeed, some galaxies of the sample can be 
considered as sosie of another. Note that each calibrating 
Cepheid has at least itself as a sosie. Obviously, we will not 
consider this special case. We will accept two Cepheids as 
sosie when the difference of their log P is smaller than 0.07. 
With a PL slope of « —3, this will give an uncertainty 
« 0.2 mag. in the distance modulus. We adopt the ratio 
Rv/Ri = 1-69 because it corresponds to the most widely 
accepted one (it corresponds to a ratio of total-to-selective 
absorption Av/{Av - Aj) = Rv/{Rv - Ri) = 2.45). 

In tableH we give the distance moduli obtained with 
equation |11| for 23 Cepheids which are sosie of another 
calibrator. In Figure 1 the comparison of the calculated 
distance moduli with the calibrating ones is given. 

From a direct regression we find that the slope is not 
different from one (1.00 ± 0.03). The observed mean dif- 
ference between the calculated distance modulus and its 
standard value is obtained together with its standard de- 
viation: 



{^sosie l^standard) 



-0.00 ±0.24. 



(12) 



The method does not introduce any systematic shift in 
the zero point. This means that the calibrating Cepheids 
constitute a coherent system (at least for the 28 Cepheids 
used in the test). The observed standard deviation (0.24) 
is in agrement with the expected standard deviation 0.2. 



Table 4. Distance moduli we obtain with equation [Tl| for 
28 Cepheids which are sosie of another calibrator. Col- 
umn 1: Name of the galactic Cepheid; Column 3: Dis- 
tance modulus calculated from Eq, 11; Column 3: Stan- 
dard deviation on the calculated distance modulus; Col- 
umn 5: Number of sosies (excluding its own case). 



Cepheid 




std.dev. 


n 


T Vel 


10.03 


0.08 


2 


BF Oph 


9.63 




1 


CV Moil 


11.07 


0.20 


2 


V Cen 


8.93 




1 


T T O 

U bgr 


8.65 




1 


BB bgr 


9.46 




1 




11.17 


0.20 


4 


V340 JNor 


11.22 


0.21 


4 


o IN or 


9.82 


U.24 


2 


TTTT T\/f,,,-. 

U U Mus 


12.59 


0.22 


3 


U Nor 


10.73 


0.40 


5 


BN Pup 


13.00 


0.23 


3 


T ^ Pun 

IjO -L up 




nil 


■3 
O 


VW Cen 


13.14 


0.19 


2 


RY Sco 


10.81 


0.13 


6 


RZ Vel 


10.98 


0.17 


6 


WZ Sgr 


11.34 


0.18 


6 


VY Car 


11.46 


0.19 


3 


WZ Car 


13.04 


0.20 


5 


VZ Pup 


13.30 


0.19 


6 


SW Vel 


11.97 


0.21 


6 


T Mon 


10.73 


0.32 


4 


RY Vel 


12.19 


0.31 


2 


AQ Pup 


12.32 


0.15 


3 


KN Cen 


13.19 


0.08 


3 


L Car 


8.65 


0.09 


2 


U Car 


11.38 


0.46 


3 


SV Vul 


11.39 




1 



4. Application to extragalactic Cepheids 

Preliminary determination of extragalactic distance 
moduli 

The method is applied to the 1840 Cepheids of Table |. To 
accept two Cepheids as sosie, we still adopt the criterion 
|logP — log Pel < 0.07 which guarantees that the stan- 
dard deviation is about 0.2 mag., assuming a PL slope of 
—3. We adopt the ratio Ry/Ri = 1.69 which corresponds 
to the first order terms proposed by Caldwell & Coulson 
(1987) and Laney & Stobie (1993). This is also the value 
adopted by Freedman et al. (2001), following Cardelli et 
al. (1987), for their HST key project about Cepheids ^ 
For each of the 36 host galaxies we plot the different dis- 
tance moduli given by Eq. Ill] as a function of log P. This 
result appears in Figure ^. 

The most important feature to point out is a signif- 
icant trend leading to higher distance moduli for long 



* This value corresponds to a ratio of total-to-selective ab- 
sorption Av/iAv - Ai) = Rv/{Rv - Ri) = 2.45. 
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Fig. 1. Comparison of standard distance moduli with 
those calculated from the method of sosie. The solid line 
corresponds to a slope of one and a zero-point of zero. 
Open circles represent the points for which there is only 
one determination and then no standard deviation. 

period Cepheids. This trend is visible for almost all the 
host galaxies. This is visible even for nearby galaxies if 
short periods are observed. For distant galaxies the trend 
is visible also at long periods. This was expected from the 
incompleteness bias we discussed elsewhere (e.g., Lanoix 
et al., 1999a). Another signature of the bias comes from 
the fact that only nearby galaxies (IC1613, IC4182, LMC, 
NGC224, NGC3109; NGC5253) have Cepheids with short 
periods. This clearly depends on the limiting magnitude 
of the considered host galaxy. This important question is 
discussed in the following subsection. 

4-2. Correction for the incompleteness bias 

In order to get the proper distance moduli we have to 
correct for the incompleteness bias. In a previous paper 
(Lanoix et al. 1999a) we suggested using a rule of thumb to 
avoid this bias. The rule consists of using only logP values 
larger than a given limit logPi- This limit is expressed as: 



logPi 



Vh 



6- 2cr 



(13) 



Unfortunatelly, this method does not take into account 
the pieces of information contained in smaller periods. The 
detailed theory of this incompleteness bias was given by 
Terrikorpi (1987) in the study of galaxy clusters. The bias 
for extragalactic Cepheids is of the same nature because 



the Cepheids of a given galaxy are all at the same distance 
from us, like the galaxies of a cluster. Assuming that the 
dispersion ct at a given logP is constant, the basic equa- 
tions adapted to the problem of extragalactic Cepheids are 
the following (for the sake of simplicity we will consider 
only the V band): 

The observed distance modulus /i will appear smaller 
than the true one. The bias A/i at a given logP is: 



A/i 



A 



TT 1 

where 



erJ{A) 



Oy log P 



aV2 



and 



dt. 



(14) 



(15) 



(16) 



In these equations a^, and are the slope and zero 
point of the PL relation. We adopt the values found by 
CFG from their galactic sample, i.e., Oy = —3.037 and 
by — —1.021 ^. Note that this requirement seems to re- 
duce the interest of the sosie method because the slope 
and zero-point are needed anyway. In fact, the slope and 
zero-point appear only as parameters in a second order 
correction. 

Two additional quantities are required to apply these 
equations: 

— The limiting magnitude Vjim- 

— The standard deviation a of the PL relation at a con- 
stant logP. 

The first quantity is derived from the histograms of {V) 
presented in Figure || for each galaxy. We adopt for Vum 
the fainter edge of the most populated class. In a few cases 
where the histogram has no dominant class, we move the 
value by ±0.5 magnitude. Vum values do not change sig- 
nificantly when one changes the binning size. Only one 
galaxy (NGC5457) changed by more than the binning size, 
but its histogram shows two classes with almost the same 
population. Nevertheless, the global influence of a change 
in Viim is discussed in section 4.3 (Table ^ and we show 
its influence on each individual galaxy in Table |^. The sec- 
ond quantity (tr) is derived by a direct linear regression 
on each plot of Figure |^. The adopted quantities Vum and 
(7 are listed in columns 2 and 3 of Table ^. 

These parameters being fixed, there is no free parame- 
ter to adjust the bias curve to the plot of Figure |^ except 
the distance modulus /i itself which is then determined 
through an iterative process. The final bias curves are plot- 
ted in Fig.|^ for each host galaxy. In column 9 of Table || 
we give the number of remaining sosies after the cut-off at 
Vurn ■ In Fig.^ the points which are rejected by the cut-off 
are represented by crosses. 

^ GFG give bv = —4.058 because they consider the zero-point 
at log P = 1 



Paturel G., et al.: Calibration of the distance scale 






ID ID 




ID ID 




ID ID H U 







i 








■ . . . 1 . . . . 1 . ... 













































—M ■ 1 1 1 1 ■ ■ 1 1 ■ ■ 1 ■ 1 l_ 









I. 





>3 II 






"H II 





Paturel G., et al.: Calibration of the distance scale 





\ ^ 



A 1 u 


" j it.' : 













Paturel G., et al.: Calibration of the distance scale 



9 



4-3. Analysis of the results 

Freedman et al. (2001) recently published their final study 
of their HST keyproject (HSTKP). They publish distance 
moduli calculated differently to those used here. They 
calibrate the PL relation with the LMC distance mod- 
ulus, assumed to be ^{LMC) = 18.5. They, adopt the 
V- and I-band PL relations and an extinction law giving 
Rv/Ri = 1.69. In order to avoid bias, they cut their sam- 
ple at a given limiting period logP; as explained above 
and they apply a small (but still uncertain) correction for 
metallicity effect. 

The comparisons between the HSTKP results and our 
solution is shown in Fig.^ for 31 galaxies in common. 
There is a fair agreement. A direct regression between 
HSTKP distance moduli and ours leads to a slope which 
is not significantly different from one (1.017 ± 0.010) and 
a zero point difference which is not significantly different 
from zero (—0.11 ± 0.16). Assuming both determinations 
carry the same uncertainty, this means that our distances 
are good within 0.16/\/2 = 0.1 magnitude. This is the 
internal uncertainty. 

From a detailed check of Figj^ one sees a slight de- 
parture from a slope of one at large distances. The effect 
is, on average, 0.17 mag. for ^ larger than 30 mag. Two 
possibilities can explain this discrepancy: 

— The PL relation of the GFG sample shows a departure 
from linearity for large log P. This effect is visible (see 
for instance Figure 4 in GFG) even when one excludes 
the three overtone Cepheids (logP < 0.6). Judging by 
the error bars of individual points, this non-linearity 
seems real. 

— The distance moduli of Freedman et al. may suffer from 
a small residual incompleteness bias. Using a simula- 
tion we have shown that it is difficult to remove the 
bias just by cutting the sample at a given log Pi. If we 
refer to our Fig. 7 in Lanoix et al. 1999a, one can see 
that at large distances (/i > 32) the bias may reach 
0.17 magnitude after the logP cutoff. At intermediate 
distances (29 < /i < 32) the bias may still reach 0.08 
mag. 

Three external sources of uncertainty come from: (i) 
the adopted ratio Ry/Ri, (h) the adopted limiting mag- 
nitude Viim and (iii) from the adopted PL relation used 
for second order bias correction. In order to check the 
stability of the solution, we repeated the previous calcula- 
tions with another PL relation (the one found by GFG for 
LMC), with a variation of Ry /Ri by ± 0.2 and a variation 
of Viim by zL0.5mag. The results are summarized in Table 
^, where we give the mean shift between distance moduli 
from different solutions and the adopted mean distance 
moduli (reference solution). One can see that the choice 
of the PL relation has no actual influence on the result. 
However, a change of Ry/Ri by ± 0.1 may change the 
mean distance modulus by nearly 0.2 magnitude and a 



Table 5. Test of the stability of the results. We give the 
departure from our reference solution for: 1) a different 
PL relation (note that this PL relation is used only for 
the second order bias correction 2) several Ry / Ri ratios. 





Rv/Ri 


a„ 


bv 


fj, — Href 


0.0 


1.69 


-2.769 


-4.063 


+0.03 ± 0.05 


0.0 


1.89 


-3.037 


-4.058 


-0.22 + 0.10 


0.0 


1.79 


-3.037 


-4.058 


-0.12 + 0.06 


0.0 


1.69 


-3.037 


-4.058 





0.0 


1.59 


-3.037 


-4.058 


+0.17 + 0.09 


0.0 


1.49 


-3.037 


-4.058 


+0.45 + 0.21 


-0.50 


1.69 


-3.037 


-4.058 


+0.20 + 0.27 


-0.25 


1.69 


-3.037 


-4.058 


+0.09 + 0.12 


+0.25 


1.69 


-3.037 


-4.058 


-0.05 + 0.08 


+0.50 


1.69 


-3.037 


-4.058 


-0.08 + 0.10 




Fig. 4. Comparison of the distance moduli from Freedman 
et al., 2001 and those from this paper. The general agree- 
ment is satisfactory but at large distances our distances 
become larger. 



change of Vum by 0.5mag may produce similar change. 
The influence of Vum depends clearly on the actual dis- 
tribution of magnitudes. For some galaxies the effect is 
negligible while it is large for some others. In order to give 
a better judgement of the stability of the distance modu- 
lus with respect to the adopted Vnm , we give the changes 
A/i~ when Vum is reduced by 0.5 mag (respectively, A/i+ 
when Viim is augmented by 0.5 mag). 

The actual uncertainty (internal plus external) can 
thus reach 0.3 magnitude and may be more if our actual 
sources of uncertainty act in the same sense. 



5. Conclusion 

The distance scale can be calibrated using galactic 
Cepheids. LMC provides us with numerous Cepheids lo- 
cated at the same distance. This gives a way to derive an 
accurate slope for the Cepheid PL relation. But its low 
metallicity (with respect to most of the galaxies of the 
sample) is a cause of suspicion; we are not sure that this 
slope can be applied to all kinds of metallicity. 
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Table 6. Distance moduli calculated from this paper using the ratio Ry/Ri = 1.69. Column 1: Name of the host 
galaxy. Column 2: The adopted limiting magnitude Vum- Columns 3: Standard deviation a. Column 4: The 
adopted distance modulus and its mean error. An asterisk marks the distance moduli of galaxies having nearly a Solar 
metallicity. Column 5: The change A/i~ of the distance modulus when Vum is reduced by 0.5 mag (i.e., a brighter 
limit). Column 6: The change A/x+ of the distance modulus when Vum is augmented by 0.5 mag (i.e., a fainter limit). 
Column 7: The number of sosie Cepheids after the Vum cut-off. 
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So, we preferred, in a first step, to calibrate the 
distance scale by using accurate distances of galactic 
Cepheids published by Gieren, Fouque and Gomez (1998). 
These distances are based on the geometrical Barnes- 
Evans method. 

Further, we applied the concept of 'sosie' (Paturel, 
1984) to extend distance determinations to extragalactic 
Cepheids without having to know either the slope or the 
zero-point of the PL relation. The distance moduli are ob- 
tained in a straightforward way. For the calibrating galac- 
tic Cepheids we checked the internal coherence from the 
same method. 



The correction for the extinction is made by using two 
bands {V and I) according to the principles introduced by 
Preedman and Madore (1990). There is no need for color 
excess estimation. 

Finally, the incompleteness bias is corrected according 

to the precepts introduced by Tcerikorpi (1987). Without 
any free parameters (except the distance modulus itself), 
the bias curve calculated for each individual host galaxy 
fits very well the observed distance moduli. This gives us 
confidence in our final distance moduli. Nevertheless, the 
small departure from the measurements published recently 
by Freedman et al. (2001) at distances larger than lOMpc 
{jjL = 30) must be clarified. 
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In order to bypass the uncertainty due to metallicity ef- 
fects it is suggested to consider only galaxies having nearly 
the same metallicity as the calibrating Cepheids (i.e. Solar 
metallicity). In Table ^ the distance moduli that can be 
considered as more secure are noted with an asterisk (*). 
Galaxies with A/i larger than « 0.3 mag. or with small 
n do not receive this flag. For a given ratio Ry/Ri, the 
uncertainty of the distances is about 0.1 magnitude but 
the total uncertainty may be about 0.3 magnitude. The 
choice of a given Ry/Ri ratio is a first source of uncer- 
tainty. The actual ratio depends on the extinction law in 
our Galaxy, on the extinction law in the host galaxy and 
on the color of the considered Cepheid. For the future it 
would be interesting to search for a clue allowing us to 
decide which value is the best in a given direction for a 
Cepheid in a given host galaxy. The proper determination 
of the limiting magnitude of the sample is a second source 
of uncertainty. It can be accurately determined only when 
a large number of Cepheids is available to provide us with 
good statistics. 

Presently, the calibration of the distance scale can 
barely be better than = 0.3 magnitude. Thus, the un- 
certainty on the Hubble constant, a{H) « cr^iJ/Sloge, 



cannot be better than about lOkr 



.Mpc 



Appendix A: The extragalactic Cepheid database 

The description of this database was given by Lanoix et 
al. (1999b). Because the database is no longer available 
on the world-wide-web the present data are published in 
electronic form in the A&A archives at CDS. All the data 
are made available, even when they are not used in the 
present paper, where only Normal Cepheids in V and I- 
bands are considered. Additional measurements were col- 
lected including the LMC ones by Udalski et al. (1999)0 
and those by Gibson et al. (1998, 1999). Data are now 
available for 2449 Cepheids of 46 galaxies (instead of 1061 
Cepheids of 33 galaxies). 

The identification of a Cepheid is given on a first line 
as follows: 

— the name of the host galaxy, 

— the name of the Cepheid and the bibliographic code 
from where this name is taken, 

— the adopted period (in log) 

— The classification of the shape of the light curve, fol- 
lowing Lanoix et al. (1999b) 

On this first line we also give the number of measurements 
attached to this Cepheid. Note that the Cepheid name for 
LMC is simply the Cepheid number from Udalski et al., 
without the field number (SC), that was not needed here 
(only three Cepheids appear with the same number in dif- 
ferent fields: 1, 16 and 19, but they are not in our list). 



We tried to keep the Cepheid name of the first discov- 
ery. This was not always done, e.g., the names given by 
Graham (1984) are referenced as Mad87 because of the 
renumbering adopted by Madore (1987). 

On the following lines, individual measurements are 
given: 

— the magnitude, 

— the type of magnitude (mean, maximum, minimum, 
average) coded according to Lanoix et al., 

— the photometric bands (B,V,R,I ...) coded according 
to Lanoix et al. (1999b), 

— the reference code. The full reference and the associ- 
ated code appears in the references. 

A sample is given below to show how the data are 
organized. 



^ Note that we kept 720 normal Cepheids among the 1182 
available with log P > 0.5. 



IC1613 
21.36 
20.79 
20.36 
20.14 
20.50 
22.03 
21.27 
21.39 

IC1613 
16.66 
18.98 
20.71 
19.85 
19.90 

IC1613 
15.47 
19.10 
17.75 
17.14 
16.62 
17.74 
20.44 
19.09 
19.09 

IC1613 
18.62 
20.10 
21.84 
20.97 
20.87 

IC1613 
21.13 
21.70 
21.46 



VI 



mea 
mea 
mea 
mea 
max 
min 
ave 
mea 



V20 



max 
min 
ave 
mea 



V22 



mea 
mea 
mea 
mea 
max 
min 
ave 
mea 



V25 



V53 



max 
min 
ave 
mea 

3 

max 
min 
mea 



Fr88a 
B 
V 
R 
I 
B 
B 
B 
B 

Fr88a 
H 
B 
B 
B 
B 

Fr88a 
H 
B 
V 
R 
I 
B 
B 
B 
B 

Fr88a 
H 
B 
B 
B 
B 

Fr88a 
B 
B 
B 



0.7480 
Fr88a 
Fr88a 
Fr88a 
Fr88a 
Sa88a 
Sa88a 
Sa88a 
Sa88a 

1.6220 
Ala84 
Sa88a 
Sa88a 
Sa88a 
Sa88a 

2.1650 
Ala84 
Fr88a 
Fr88a 
Fr88a 
Fr88a 
Sa88a 
Sa88a 
Sa88a 
Sa88a 

0.9600 
Ala84 
Sa88a 
Sa88a 
Sa88a 
Sa88a 

0.5900 
Car90 
Car90 
Car90 



B+ 
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